Post-merger electromagnetic emissions from disks perturbed by binary black holes 
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We simulate the possible emission from a disk perturbed by a recoiling super-massive black hole. 
To this end, we study radiation transfer from the system incorporating bremsstrahlung emission 
from a Maxwellian plasma and absorption given by Kramer's opacity law modified to incorporate 
blackbody effects. We employ this model in the radiation transfer integration to compute the 
luminosity at several frequencies, and compare with previous bremsstrahlung luminosity estimations 
from a transparent limit (in which the emissivity is integrated over the computational domain and 
over all frequencies) and with a simple thermal emission model. We find close agreement between 
the radiation transfer results and the estimated bremsstrahlung luminosity from previous work for 
electromagnetic signals above 10 14 Hz. For lower frequencies, we find a self-eclipsing behavior in 
the disk, resulting in a strong intensity variability connected to the orbital period of the disk. 

PACS numbers: 



I. INTRODUCTION: 

The near possibility of detecting gravitational waves to 
study astrophysical systems has spurred efforts to under- 
stand such systems theoretically to aid in the detection 
and interpretation of obtained signals (see, for instance, 
P, [3, H, 0, IE @] and references cited therein for further 
details and examples). Systems that radiate strongly in 
both electromagnetic and gravitational wave bands are 
particularly interesting since the combined information 
will provide unprecedented access to a number of rich 
phenomena (e.g. @, H, Q-) Electromagnetic radiation 
from compact objects can be (partially or completely) 
absorbed and scattered by dust and other surrounding 
material. Gravitational radiation, on the other hand, is 
generated in the central engine, and these waves do not 
interact with dust or material in the surrounding envi- 
ronment. 

One such system involves the collision of a binary black 
hole within a circumbinary disk [lEEILEll- For a generic 
black hole binary, the gravitational waves are radiated 
asymmetrically. As these waves carry energy and mo- 
mentum, the resultant black hole can have a kick or re- 
coil velocity [H, El EE EE E3, El , which in turn will 
induce shocks in the disk. A second source of pertur- 
bations in the disk comes from the loss of gravitational 
mass with the emission of gravitational radiation, weak- 
ening the gravitational pot ential felt by the disk. These 
effects, as discussed in [19], |2(1 [2l|, can induce possible 
emissions and some of these have been investigated in 
part in [H HI, IHHE El- These works stud y thc dy- 



namics of the disks when either mass reduction (due to 
the energy lost by the system) or recoil velocity (due to 
the asymmetric flux of gravitational radiation) effects are 
included and present preliminary estimates of emissions 
from the system with the goal of understanding possible 
observations in the electromagnetic spectrum. 



In the present work, we further investigate electromag- 
netic radiation from shocked disks by directly comput- 
ing the electromagnetic emission using radiation trans- 
fer. We build on our earlier work on disks [23|, hereafter 
referred to as Paper I, where we followed the evolution of 
disks when thc central black hole was given different re- 
coil velocities. We use this data for the evolution of disks 
to calculate the electromagnetic spectrum using radia- 
tion transfer. In this method, we assume simple models 
for the emission and absorption of electromagnetic radia- 
tion and solve for the intensities by integrating the radi- 
ation transfer equation along geodesies of the spacetime. 
We examine multiple radiation models to investigate how 
different choices in the emission models affect the time- 
dependent luminosity. One model assumes only classical 
thermal emission, and a second model incorporates both 
bremsstrahlung and thermal emission according to op- 
tical depth. Finally, wc also compare our results with 
the simple bremsstrahlung estimate that we obtained in 
Paper I. 
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II. RADIATION TRANSFER 

Electromagnetic radiation follows geodesies of the 
spacetime. As the photon is massless, the radiation fol- 
lows null geodesies, meaning that the proper separation 
between all points on the geodesic is zero. We calculate 
the geodesies from the spacetime metric g a p, and use the 
convention that Greek indices take the values {0, 1, 2, 3}. 
The null geodesies are found by solving the first order 
system 



dx>* 
~dX 

dA 



(1) 
(2) 



where r Q M7 are the Christoffel symbols, and and p M 
are the coordinates and momentum of the photon, re- 
spectively. Since we are dealing with null geodesies we 
parameterize the solution in terms of the affine parameter 
A. We solve the geodesic equations using the standard 
Runge-Kutta fourth order integration method, numeri- 
cally integrating the geodesies backwards in time, per- 
pendicularly from a viewing window. This window is set 
at different locations, in order to provide different view- 
ing angles of the disk. 

After obtaining the geodesic paths in the spacetime, we 
can now solve for the electromagnetic radiation intensity 
along these paths. The radiation transfer equation for 
the photon intensity at frequency v is (25j 



dA 



-p a u a [rjo - xoL] 



(3) 



where xo and ?7o are, respectively, the absorption and 
emission coefficients in the rest frame with respect to the 
fluid, and u a is the four-velocity of the fluid. In vac- 
uum regions, the fluid four- velocity corresponds to the 
four-velocity of coordinate observers. This first order ra- 
diation transfer equation has the solution p5| 



,(A) = /,(Ao) e / ^ Xo(A>o)UapQdA ' - 



f e J$ Xo(X'\,o)u^dX" m{x ^ VQ)UardX > (4) 

where Uq is the photon frequency at emission, and Ao 
indicates the affine parameter value where the geodesic 
begins. The integration constant I„(Ao) takes into ac- 
count the photon intensity where the geodesic begins. 
In this work, the geodesies begin far enough behind the 
disk so that there is no matter to act as a source for 
emission, and we set I„(Ao) to be zero. The fluid vari- 
ables come from a numerical solution that is specified on 
a finite difference grid. When integrating Eq. ([?]), these 
fluid variables arc interpolated both in space and time. 



As mentioned above, we employ two different emis- 
sion models in calculating the spectrum from the kicked 
disk. The first model assumes that emission from the 
optically thin region is thermal bremsstrahlung, and the 
optically thick region is a classical blackbody [26|. The 
second model uses only thermal emission. Both models 
use Kramer's opacity law for the absorption. We expect 
the first model to be more realistic, as the optically thin 
regions of the disk are expected to radiate primarily with 
bremsstrahlung, rather than thermal emission, and the 
second model provides a useful comparison to gauge the 
importance of the bremsstrahlung emission. We refer to 
these models as bremsstrahlung-blackbody and thermal 
emission models, respectively. 

The bremsstrahlung-blackbody radiation transfer 
model uses the standard bremsstrahlung cmissivity and 
a modified form of Kramer's opacity law for the absorp- 
tion. The bremsstrahlung emissivity is modeled on a 
Maxwellian plasma [27j 



r)o = 6.8 x 1CT 38 Z 2 n e n t T- 112 
erg 
s cm 3 Hz ' 



xg s (u,T)e- x 



where 



x = hv jkT. 



(5) 



(6) 



In these equations h is Planck's constant, k is Boltz- 
mann's constant, T is the temperature of the fluid, n e 
is the electron number density, rit is the ion number den- 
sity, Z — rie/rii, and gg is the temperature averaged 
Gaunt factor for free-free transitions in pure hydrogen. 
The electron and ion number densities are given by 



(7) 



where 



fi e m p 



1 + X' 



fJ-iTHp 



1 + 3X 



(8) 



Here p is the rest density of the fluid, m p is the mass of 
the proton, and X is the relative abundance of hydro- 
gen in the universe, which we set to be X = 0.75. We 
use a temperature averaged Gaunt factor, gg, based on 
[27| . where different approximations are given depend- 
ing on the values of T and v. To completely specify gs, 
some criteria must be used to set the transition between 
approximations valid for v -C kT/h and v sa kT/h in 
the regime T(K) > 3 x 10 5 Z 2 . We set the transition at 
the values \v, T] for which gg in the two approximations 
coincide, giving 
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and T < 3 x 10 5 Z 2 



and T > 3 x 10 5 Z 2 



(9) 



where T is in units of Kelvin, and v in Hertz. Follow- 
ing [26| . we use for the absorption a modified Kramer's 
opacity law (compare to Eq. [12] below) , 



Xo = 5 x io 2 y r 



-7/2 



1 



(10) 



This modification has the net effect of adding thermal 
radiation for optically thick regions, while optically thin 
regions radiate with bremsstrahhmg. In optically thick 
regions, all photons arc absorbed and dI v /dX — > 0, giving 
the expected blackbody intensity 



I v = — ocTV 5 - 

Xo l-e~ 



(11) 



In the second radiation transfer model we consider only 
the emission of thermal radiation. We use Kramer's opac- 
ity law for the absorption 

Xo = 5 x 10 24 p 2 T- 7/2 cm-\ (12) 

and specify the emissivity using Kirchoff 's law [28| 



r)o = XoB v (T), 
where B V (T) is given by the Planck law 

2/w 3 /c 2 



Bu(T) 



1 



(13) 



(14) 



Thermal radiation is emitted by matter in thermody- 
namic equilibrium, while blackbody radiation occurs 
when the radiation is itself in thermodynamic equilib- 
rium. In optically thick regions (dl^/dX — > 0), wc again 
reduce to the blackbody intensity, I® B = B V {T). 

Having specified our emission models, we briefly men- 
tion some physical effects not yet contemplated in our 
analysis. As photons move away from the gravitational 
sources, their frequencies change due to gravitational red- 
shift. As the redshift varies depending on the location at 
which the photon is emitted, this leads to a broadening of 



the emission lines. A simple estimate in a Schwarzschild 
spacetime shows that gravitational redshift accounts for 
typical frequency changes of 5% percent or less if we con- 
sider photons being emitted from a minimum distance to 
the black holes of about 10 Rs, consistent with the disk 
in our model. Broadening will occur also due to relativis- 
tic Doppler, especially in the inclined views. Considering 
the fluid velocities in our model, up to about 0.28c, we 
estimate a maximum Doppler shift of about 25%. As 
the frequency bins used here are much larger than those 
changes, we do not include the line broadening caused 
by Doppler or gravitational redshift in our model, and 
we work only with the original emission frequencies. 

Finally, in Paper I we estimated the total 
bremsstrahlung luminosity by assuming a transpar- 
ent gas. In this approximation the luminosity is 
obtained by integrating the bremsstrahhmg emissivity, 
tB cx P 2 T 1 ^ 2 , over the computational domain. Wc 
compare the results computed here with this earlier 
model, and find that this simple estimation works very 
well for the higher frequencies that we consider. 



III. RESULTS 

In paper I the dynamics of a disk affected by a possible 
recoil velocity of a newly merged black hole was studied. 
A salient feature of the analysis presented there is that 
the recoil influences strongly the disk's behavior as long 
as it has a component on the orbital plane. This induces 
shocks in the disk which in turn heat the gas and can 
lead to emissions with a variability tied to the orbital 
period at the disk. The results presented in Paper I were 
obtained assuming the simple ideal fluid equation of state 



P = (r - l)pe, 



(15) 



where P is the fluid pressure, e is the specific internal en- 
ergy, and we choose T = 5/3, considering the gas as being 
monatomic. Wc assume that there is no energy trans- 
fer from the fluid to the electromagnetic radiation, and 
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this equation of state allows for fluid flows, e.g. shocks, 
that are not isentropic. To calculate the emissivity and 
absorption coefficients, we obtain the fluid temperature 
from the ideal fluid law as 



T = 



fim p 



where 



f 

/' 



1 

Me 



1 



(16) 



(17) 



and p e and pn are given in Eqs. ([8]). 

The simulations in Paper I correspond to a fixed space- 
time background, and thus some variables can be inde- 
pendently rescalcd a posteriori. For example, one can 
rescale the black hole's mass and the disk's density in- 
dependently, which in turn fixes other quantities such 
as the luminosity. This scale invariance, however, does 
not hold for the emission and absorption coefficients in 
the radiation transfer equation ([4J , and we set the phys- 
ical scales in the following manner. We start by setting 
a central black hole mass Mbh = 1 x 10 s Mq, consis- 
tent with super-massive black holes at galactic centers. 
This choice fixes the location of the maximum fluid den- 
sity p m in our model to r m sa 4.5 x 10 9 km, approxi- 
mately 15 Schwarzschild radii (Rs) from the center of the 
disk. To estimate a physically reasonable value for the 
disk's density, we use an alpha-disk model [2!| . Note that 
this model, unlike ours, corresponds to an accreting disk. 
However, as an estimate we assume that the density in 
the alpha-disk model will be similar to our disk at the pre- 
merger stage, before the inner edge of the disk "freezes" 
and accretion stops 0, EH • We set the accretion rate in 
the alpha disk model to be an Eddington-limited accre- 
tion rate of 1 M & /yea.r. We choose a = 0.1, a value typ- 
ical for cataclysmic variables, low-mass X-ray binaries, 
and black hole transients. Using these parameters in the 
alpha-disk model we obtain a density p m 0.01 g/cm 3 
at r = r m . Hence, we scale the disk density so that 
p m = 0.01 g/cm 3 . 

As discussed in Paper I, the qualitative behavior for 
all of the disks that we considered was similar, provided 
that the recoil velocity has a component orthogonal to 
the disk's angular momentum. We thus adopt as a rep- 
resentative example the data obtained for a kick perpen- 
dicular to the axis of rotation at 3000 km/s. We place a 
square-shaped window at three different locations, pro- 
viding viewing angles of 0°, 60°, and 75° with respect to 
the disk's rotational axis, while the azimuthal direction is 
chosen perpendicular to the kick. In all cases the window 
is located at a distance of 1.05 x 10 10 km (~ 35i?s) away 
from the black hole (at the edge of the computational 
domain employed in Paper I). The area of the observa- 
tion window was 4.41 x 10 20 km 2 so that the whole of 
the disk could be analyzed in such window. The num- 
ber of geodesies used for the window and the number of 
points used to integrate the radiation transfer equation 



were selected so that the results did not show appreciable 
differences upon further increasing them. 

We concentrate on photons having initial frequencies 
from infrared to gamma ray (within the range 10 12 - 
10 21 Hz) and integrate the radiation transfer equation to 
the window where we examine the resulting luminosity vs 
initial frequency. We compare results obtained from the 
bremsstrahlung-blackbody and thermal models using ra- 
diation transfer. We also compare brcmsstrahlung lumi- 
nosity estimated in the transparent limit (Paper I) with 
the results from radiation transfer. This allows us to ex- 
amine interesting frequency dependencies not highlighted 
by the brcmsstrahlung luminosity estimate by itself. 

For instance, Figures [T] and [5] illustrate the lumi- 
nosity for the case of soft x-ray photons (initial fre- 
quency of 10 18 Hz). Figure Q] shows normalized lu- 
minosity curves (with respect to the initial value) ob- 
tained via the radiation transfer equation for a win- 
dow located at 0° viewing angle (the top view) while 
Figure [2] shows normalized luminosity curves for dif- 
ferent viewing angles. In Figure [T] the results for two 
radiation models are shown: bremsstrahlung for emis- 
sivity and modified Kramer's opacity law for absorp- 
tion (bremsstrahlung-blackbody), and thermal emissivity 
with unmodified Kramer's opacity law (thermal). Also 
included is the normalized bremsstrahlung luminosity 
from Paper I which was calculated assuming a transpar- 
ent gas, and integrated over all frequencies (see Eq. 16 
of Paper I). The good agreement seen in Figure [T] be- 
tween the bremsstrahlung-blackbody radiation transfer 
luminosity curve and the bremsstrahlung luminosity es- 
timate from Paper I indicates that the disk is largely 
transparent at 10 18 Hz. 

Interestingly this behavior no longer holds for lower 
frequencies. For instance, Figure [3] illustrates the results 
obtained in the infrared (initial frequency of 10 12 Hz). At 
this lower frequency, much of the disk is opaque and a 
self-eclipsing behavior is increasingly evident as the incli- 
nation angle increases. More importantly, after the tran- 
sient behavior, a strong increase in relative luminosity 
is observed as shocks significantly affect the dynamics of 
the disk. Finally, Fig. Q] shows the electromagnetic spec- 
trum calculated from 36 frequencies in the range from 
10 12 -5 x 10 21 Hz. 

As an illustration of the disk structure observed at dif- 
ferent initial frequencies and viewing angles, Fig. [5] and 
Fig. [HI illustrate the intensity obtained at different times. 
The color map shows areas of higher intensity emission. 
Comparing these two figures one can clearly see the dif- 
ferences between a mostly transparent case (Fig. [5]) , and 
one with absorption (Fig. [5]), showing self eclipsing, which 
supports the oscillations shown in Fig. [31 seen only in the 
opaque case. 
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FIG. 1: Normalized luminosity as a function of time using 
different emission models for a frequency of 10 18 Hz. The hor- 
izontal axis shows time in seconds after the kick. In all cases, 
the results are normalized with respect to their initial value. 
The curve labeled "bremsstrahlung estimate" (dash-dot line) 
is calculated by assuming a transparent gas and integrating 
p2yi/2 over computational domain. The curve labeled 
"thermal" (dashed line) is obtained using the thermal radia- 
tion model described in the text for a 0° viewing angle (top 
view). The one labeled "bremsstrahlung-blackbody" (solid 
line) is obtained by using the bremsstrahlung-blackbody ra- 
diation model described in the text, also at a 0° viewing angle. 
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FIG. 2: Normalized bremsstrahlung-blackbody luminosity as 
a function of time for different viewing angles for a frequency 
of 10 18 Hz. The horizontal axis is time in seconds after the 
kick. 



IV. FINAL COMMENTS 

We have analyzed the possible emission of a disk per- 
turbed by a recoiling black hole by directly imaging 
the disk with radiation transfer. For frequencies above 
~ 10 13 Hz our studies reveal a characteristic luminos- 
ity variability that could be present in observations of 



FIG. 3: Normalized luminosity as a function of time for differ- 
ent viewing angles for a frequency of 10 12 Hz. The luminosity 
was normalized with respect to its initial value. The hori- 
zontal axis labels time in seconds after the kick. At 10 12 Hz 
much of the disk is opaque and exhibits self-eclipsing behav- 
ior. The oscillations seen after t = 1 X 10 6 seconds are tied 
to the orbital period of the disk and are more pronounced as 
the viewing angle increases. 



this kind of systems. This variability shows the same 
main features as those in our previous results, in which a 
simpler model -without radiation transfer- was used to 
estimate the luminosity. 

For frequencies below ~ 10 13 Hz the spectrum matches 
the Rayleigh- Jeans frequency dependence expected from 
a blackbody. Furthermore, at these lower frequencies, a 
self-eclipsing behavior at lower frequencies is clearly seen 
as the disk orbits around the black hole, so that bright 
and opaque regions alternatively superpose along the line 
of sight, resulting in an induced variability tightly tied 
to the orbital motion and inclination. This effect is less 
pronounced as higher frequencies are considered and so 
it is unlikely this feature will be observed. 

Nevertheless, a strong overall variability is observed 
at all frequencies associated with shock heating in sec- 
tors of the disk which could be detected provided that 
the increase in the disk's luminosity prior to the recoil 
is sufficiently low. This appears as a likely scenario as 
the binary hollows out the inner regions of the disk [T^j • 
One should be cautious that uncertainties on the details 
of these disks and emission processes can affect the char- 
acteristics of the luminosities obtained. In particular, 
the work of (24|, which employs a different equation of 
state, displays interesting differences with the ones pre- 
sented here. However, regardless of particular differences 
induced by details of the equation of state employed, a 
common message of all related works [2(J Hi], [22|, HH, [24| 
is that emissions from circumbinary disks around merg- 
ing black holes possibly display particular features that 
could be detected. These features would follow the black 
hole merger in the time-frame of months and would re- 
quire a good localization in the sky to follow gravitational 
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wave observations with these electromagnetic counter- 
parts. As already discussed in the literature (e.g. [U and 
references cited therein) LISA's localization will unlikely 
be able to single out a preferred galaxy and so prospects 
of detection will be greatly enhanced if pre-cursors or 
prompt emissions in the electromagnetic wave band are 
also available 0] . In the case of binary black hole merg- 
ers that give rise to the scenario analyzed in this work, 
such a p ossibility is just beginning to be analyzed in de- 
tail [M [Mill [MS]. 



FIG. 4: The electromagnetic spectrum computed 1.25 x 10 6 
seconds after the kick using the "bremsstrahlung-blackbody" 
radiation model at a 0° viewing angle. The horizontal axis is 
the log 10 of the frequency. The vertical axis is the log 10 of 
the intensity integrated over the window area. In the lower 
frequencies of the spectrum 10 12 — 10 13 Hz we see the ex- 
pected Ray leigh- Jeans v 2 frequency dependence. The small 
kink in the spectrum in between 10 18 and 10 19 Hz is due to 
the behavior of the Gaunt factor in that region. The spec- 
trum was calculated using 36 frequencies, spaced evenly in log 
space from 10 12 -5.6 x 10 21 Hz. 
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FIG. 5: The intensity for 10 18 Hz shown, from top to bottom, at t = 1.75 x 10 5 , 1.25 x 10 6 , 1.5 x 10 6 , and 1.75 x 10 6 seconds 
after the kick using the "bremsstrahlung-blackbody" radiation model. The colormap is a log scale; The intensity for all the 
images has been normalized using the maximum intensity found at t = 1.25 x 10 6 seconds after the kick. The images on the 
left /center /right columns correspond to viewing angles of 0°/60°/75° respectively. The number of geodesies used was 40,000. 
The axis scale of the images is in units of mass of the black hole. 
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FIG. 6: The intensity for 10 12 Hz at t = 1.75 x 10 s , 1.25 x 10 6 , 1.5 x 10 6 , and 1.75 x 10 6 seconds after the kick using the 
"bremsstrahlung-blackbody" radiation model. The colormap is a log scale; The intensity for all the images has been normalized 
using the maximum intensity found at t = 1.25 x 10 6 seconds after the kick. The images on the left/center/right columns 
correspond to viewing angles of 0°/60°/75° respectively. The number of geodesies used was 40,000. The axis scale of the images 
is in units of mass of the black hole. 



